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The  ground  shock  calculations  codes  currently  used  to  predict  the  states 
of  stress  and  ground  motions  induced  in  natural  earth  masses  by  explosive  deto¬ 
nations  are  deterministic  tools;  that  is,  their  input  parameters  (soil  stress- 
strain  and  strength  properties,  soil  density,  applied  airblast  surface  over¬ 
pressure,  etc.)  are  specified  as  single-valued  quantities.  In  actuality,  how¬ 
ever,  both  the  properties  of  earth  materials  and  the  characteristics  of  the 
airhlast  pulse  are  random  variables.  Consequently,  the  resulting  states  of 
stress  and  ground  motion  in  a  natural  earth  mass  (Continued) 
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20.  ABSTRACT  (Continued) 

are  also  random  variables.  Therefore,  the.  ground  shock  calculation  problem 
should  be  treated  probabilistically. 

-This  paper  describes  the  development  of  a  probabilistic  solution  for 
stress  wave  propagation  in  homogeneous  bilinear  hysteretic  materials  subjected 
to  an  exponentially  decaying  surface  airblast  pulse.  The  partial  derivative 
method  is  used  to  calculate  the  time  histories  of  the  moments  of  the  dependent 
(output)  variables,  consisting  of  particle  displacement,  particle  velocity,  and 
stress,  in  terms  of  functions  of  moments  of  the  independent  (input)  variables. 
The  partial  derivatives  are  evaluated  numerically  using  finite-difference 
approximations.  Several  sensitivity  analyses  are  conducted  to  illustrate  how 
the  solution  technique  can  be  used  to  quantify  and  rank  the  relative  effects  of 
input  variabilities  or  uncertainties  on  the  dispersion  of  output  quantities. 
Finally,  application  of  the  methodology  for  conducting  probabilistic  analyses 
of  airblast-induced  vertical  ground  shock  is  demonstrated,. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used 

units  as  follows: 

Multiply 

in  this  report  can 

By 

be  converted  to  SI  (metric) 

To  Obtain 

feet 

0.3048 

metres 

feet  per  second 

0.3048 

metres  per  second 

inches 

0.0254 

metres 

kips  (1000  lb  force)  per  square 

inch  6.894757 

megapascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 
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PROBABILISTIC  SOLUTION  FOR  ONE-DIMENSIONAL  PLANE  WAVE  PROPAGATION  IN 
HOMOGENEOUS  BILINEAR  HYSTERETIC  MATERIALS 

PART  I:  INTRODUCTION 

1.  The  ground  shock  calculation  techniques  currently  used  to  predict 
the  states  of  stress  and  ground  motions  induced  in  natural  earth  masses  by 
explosive  detonations  are  deterministic  tools — that  is,  their  input  parame¬ 
ters  (media  stress-strain  and  strength  properties,  density,  surface  airblast 
loading,  etc.)  are  specified  as  deterministic  quantities  or  functions.  In 
actuality,  however,  both  the  properties  of  earth  materials  and  the  charac¬ 
teristics  of  airblast  pulses  are  random  variables.  Consequently,  the  random¬ 
ness  of  these  input  variables  indicates  that  the  resulting  states  of  stress 
and  ground  motion  in  natural  earth  masses  are  also  random  variables. 

Therefore  the  ground  shock  calculation  problem  should  be  treated  probabilis¬ 
tically. 

2.  In  general,  airblast-induced  ground  motion  from  a  surface  burst  can 
be  analyzed  as  an  axisymmetric  two-dimensional  (2D)  problem.  And,  if  the 
ratio  of  the  propagation  velocity  of  the  airblast  pulse  traversing  the 
ground  surface  is  much  greater  than  the  propagation  velocity  of  the  stress 
wave  in  the  medium  (superseismic  conditions),  the  near-surface  ground  motions 
outside  the  crater  are  predominantly  vertical,  in  which  case  one-dimensional 
(ID)  plane  wave  calculations  are  usually  adequate  for  predicting  free-field 
response.  It  is  appropriate,  therefore,  as  well  as  useful  to  commence  an 
examination  of  probabilistic  ground  shock  in  terms  of  the  ID  problem  before 
examining  the  more  cumbersome  2D  problem. 

3.  A  simple,  yet  powerful,  deterministic  model  for  predicting  ID  wave 
propagation  phenomena  in  earth  media  is  the  analytic  solution  for  a  bilinear 
hysteretic  material  loaded  by  a  decaying  surface  airblast  pulse  developed  by 
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Salvador i ,  et  al . ,  in  1960  (Reference  1).  Although  the  real  stress-strain 
properties  of  soil  are  only  approximated  by  this  bilinear  relationship,  the 
solution  accounts  for  the  major  observed  features  of  stress  wave  propagation 
through  earth  materials — namely,  the  attenuation  of  stress  and  particle 
motion  with  depth.  This  paper  is  concerned  with  the  conversion  of  this 
deterministic  solution  into  a  probabilistic  solution  using  the  method  of 
partial  derivatives  (Reference  2) .  The  method  of  partial  derivatives  greatly 
alleviates  the  burden  of  conducting  brute  force  Monte  Carlo  analyses  for 
probabilistic  wave  propagation  problems.  Furthermore,  the  partial  derivative 
method  has  the  advantage  that  it  can  be  used  to  quantitatively  rank  the 
relative  effects  of  input  variabilities  (uncertainties)  on  the  dispersion  of 
the  output  quantities,  i.e.,  particle  displacement,  particle  velocity,  and 
stress . 
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PART  II:  DETERMINISTIC  SOLUTION  FOR 


ONE-DIMENSIONAL  WAVE  PROPAGATION  IN  BILINEAR  HYSTERETIC  MATERIALS 


Description  of  Bilinear  Model 

4.  The  bilinear  hysteretic  model  (Reference  1)  was  first  used  to 

approximate  the  stress-strain  behavior  of  soils  in  states  of  uniaxial  strain. 
The  actual  soil  stress-strain  curve  during  virgin  loading  is  approximated  as 
a  straight  line  which  is  defined  by  modulus  .*  During  unloading  and 

subsequent  reloading  to  a  previous  maximum  stress  level,  the  actual  stress- 
strain  relation  is  different  from  the  loading  relation  and  is  approximated  by 
another  straight  line,  defined  by  modulus  .  When  =  M^  ,  the  model 

corresponds  to  a  linear  elastic  material.  The  propagation  velocity  of  a 
virgin  loading  stress  wave  is  given  by 

C  =  *¥7p~  (1) 

o  o 

where  p  is  the  mass  density  of  the  material  and  the  propagation  velocity  of 
an  unloading  or  reloading  stress  wave  is  given  by 

C1  =  /M^p  (2) 

The  bilinear  hysteretic  material  is  therefore  completely  described  by  three 

material  constants:  M  ,  M,  ,  and  p  . 

o  1 

Boundary  Loading 

5.  The  dynamic  boundary  load  considered  in  Reference  1  allows  the 
solution  to  be  used  for  blast-type  problems.  It  is  a  pulse  characterized  by 

*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  A) . 
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an  instantaneous  rise  to  peak  pressure  (a  shock  front)  followed  by  an 
exponential  decay.  Its  expression  is 

P(t)  =  Pq  exp(-t/i)  (3) 

where  Pq  is  the  peak  applied  pressure,  t  is  time,  and  t  is  the  expo¬ 
nential  time  constant  (time  at  which  pressure  has  decayed  to  0.3678  P  ). 

The  airblast  pulse  is  therefore  completely  characterized  by  two  parameters: 

P  and  x  . 
o 

Response  of  the  Medium  to  Surface  Boundary  Loading 

6.  The  ID  problem  treated  by  this  deterministic  solution  may  be  viewed 
as  a  semi-infinite  body  of  bilinear  hysteretic  material  uniformly  loaded  at 
its  free  surface  by  the  pressure  pulse  P(t)  described  above.  Surface 
motion  and  stress  are  assumed  to  be  zero  before  application  of  the  load. 
According  to  Reference  1,  the  response  of  the  medium  in  terms  of  time 
histories  of  stress  a  and  particle  velocity  U  .n  the  direction  of 
propagating  wave  at  a  generic  depth  Z  is  given  as 


o(Z,t) 


U(Z,t) 
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Integration  of  Equation  5  results  in  the  following  expression  for  particle 
displacement 


u(z’t)  =  ?°{A  -  A)(c  H'fe  ■  A/]  “  “pKc  •  k)l]  1 
+ £  |2a"(t  ■  ¥ +  exp[“"(t  “  A)l]+  “p[-"(t + k)l] 

P(‘  ■  A)l]  -  •w[-n( +  k)f\ 


-  exp  I  -i 


The  parameter  a  in  Equations  A  through  6  is  given  by 


/M,/M  -  1 

1  O 

a  =  ■  ~  - 

/m./m  +  1 

1  O 


(7) 


To  account  for  the  travel  time  of  the  wave  to  the  depth  of  interest  the 
initial  value  of  t  in  Equations  4  through  6  must  correspond  to  the  arrival 
time  tfl  at  that  depth  (i.e.,  these  equations  hold  for  t  t  ) ,  which  is 
given  by 


t  =  Z/C 
a  o 


(8) 


Equations  4,  5,  6,  and  8  (which  relate  the  four  dependent  variables  a  , 

U  ,  U  ,  and  t  and  the  five  independent  variables  Mq  ,  ,  p  ,  P  ,  and 

t)  provide  the  complete  deterministic  solution  for  ID  stress  wave  propagation 

in  bilinear  hysteretic  materials.  Note  that  since  the  stress  and  particle 

velocity  maxima  of  this  solution  always  occur  at  the  wave  front,  t  in 

Equations  4  and  5  can  be  replaced  by  t  to  obtain  expressions  for  the  peak 

stress  o  and  peak  particle  velocity  U  as  a  function  of  depth, 

max  max  v 
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PART  III:  PROBABILISTIC  ANALYSIS 


7.  The  purpose  of  a  probabilistic  analysis  is  to  develop  a  method  by 
which  the  variability  or  uncertainties  in  the  independent  parameters  in  a 
particular  problem  can  be  evaluated  or  estimated  in  terms  of  their  effects  on 
the  dispersion  of  the  dependent  (output)  variables.  A  highly  useful  proce¬ 
dure  for  implementing  such  an  analysis  is  to  apply  the  method  of  partial 
derivatives,  described  in  References  2  and  3,  to  a  deterministic  solution  of 
the  problem.  The  method  gives  approximations  for  the  moments  of  the  dependent 
variables  in  terms  of  functions  of  moments  of  the  independent  variables. 

That  is  if  a  random  variable  Y  is  functionally  related  to  the  random 
variables  X. 

l 


Y  =  Y(XrX2,...,Xn) 


(9) 


and  if  the  Xi  are  uncorrelated,  then  according  to  the  partial  derivative 
method  the  approximations  for  the  expected  value  of  Y,  E[y]  and  the  vari¬ 
ance  of  Y,  Var[Y]  ,  are 


E[Y]  ~  Y(p1,p2, 


'un>  +  \ 


n 

l 


32Y 


i=l  3Xj 


Varfx^  (10) 


(p1.u2» ' 


’V 


Var [v] 


n 

l 

i=l 


VarCxJ 


(ID 


where  (p^,p2> . . . .y^)  denotes  the  respective  mean  of  the  random  variables 
(X^,X2, . . . ,Xn)  .  The  first  term  on  the  right-hand  side  of  Equation  10 
corresponds  to  the  mean  value  of  Y  ,  i.e.,  the  value  of  Y  obtained  using 
the  mean  values  of  all  of  the  random  variables.  The  second  term  represents 
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the  contributions  to  the  expected  value  of  Y  due  to  uncertainties  in  the 


X^  .  The  second  term  is  negligible  if  Var[x^]  and  the  nonlinearity  in 
the  function  Y  itself  are  not  large.  As  pointed  out  in  Reference  2, 
Equation  11  "may  be  interpreted  as  meaning  that  each  of  the  n  random 
variables  X^^  contributes  to  the  dispersion  of  Y  in  a  manner  proportional 
to  its  own  variance  Var[x^]  and  proportional  to  a  factor 


which  is  related  to  the  sensitivity  of  changes  in 


Y  to  changes  in  X^."  This  interpretation  can  be  used  to  conduct  sensi¬ 
tivity  analyses  to  quantify  and  rank  the  relative  effects  of  the  input 
variabilities  or  uncertainties  on  the  dispersion  of  the  output  quantities. 

8.  The  partial  derivatives  in  Equations  10  and  11  can  be  evaluated 
analytically  if  an  explicit  expression  is  available  for  the  dependent 
variable  Y  .  However,  as  pointed  out  by  Mlakar  (Reference  4),  in  many 
cases,  even  when  an  explicit  relation  for  Y  does  exist,  it  is  often  more 
convenient  to  evaluate  these  partial  derivatives  numerically  using  finite- 
difference  approximations.  Following  the  method  proposed  by  Mlakar 
(Reference  4),  the  partial  derivatives  may  be  expressed  as 


3Y 

3X. 

l 


(12) 


. 2V  Y (y  , . . . ,  y  -kS  ,  . . 
oi  i  1  i 


. , y  ) —  2Y (y . , .. .  ,  ,  ■ 

nil 


■  *  b^) "^"Y  (y^ . .  ^~FkS^  t  •  •  •  b^) 


i  —  1,2, ...,n 


The  first  partial  derivative  (Equation  12)  is  calculated  from  the  functional 
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values  of  Y  at  k  standard  deviation  kS^  above  and  below  the  mean 
value  of  X^  (or  p^)  where  k  is  a  constant  and  is  the  standard 

deviation  of  X^  .  Similarly,  the  second  partial  derivative  (Equation  13) 
is  calculated  from  the  mean  value  of  Y  and  the  functional  values  of  Y  at 
kS^  above  and  below  .  The  term  kS^  therefore  may  be  viewed  as  a 
finite-difference  mesh  spacing.  The  value  of  k  is  heuristically  taken  to 
be  1.0  in  Reference  4.  However,  if  either  the  standard  deviation  or 

the  nonlinearity  in  the  function  Y  is  large,  smaller  values  of  k  should 
be  tried  until  the  results  are  no  longer  affected  by  further  reductions  in 
k  .  The  probabilistic  calculations  reported  in  this  paper  (next  section) 
were  conducted  using  k  =  0.5  .  In  order  to  determine  the  sensitivity  of  the 
results  to  variations  in  k  several  calculations  were  conducted  for  k 
varying  from  1.0  to  0.1.  These  calculations  indicated  that  an  order  of 
magnitude  reduction  in  the  value  of  k  decreased  the  numerical  results  on 
the  average  by  6  percent  and  at  most  by  12  percent.  Reducing  the  value  of  k 
from  0.5  to  0.1,  however,  changed  the  results  of  the  calculations  on  the 
average  by  1  percent  and  at  most  by  3  percent.  For  most  practical  problems 
it  would  therefore  appear  that  a  value  of  k  within  the  range  of  0.1  to  0.5 
would  be  satisfactory. 

9.  For  each  of  the  dependent  variables  o  ,  U  ,  U  ,  and  tfl 

(Equations  4,  5,  6,  and  8),  Equations  12  and  13  can  be  used  to  calculate 

their  first  and  second  partial  derivatives  with  respect  to  the  five  independent 

variables  M  ,  M,  ,  p  ,  ?  ,  and  t  .  These  partial  derivatives  can  then 

o  l  o 

be  used  in  Equations  10  and  11  to  calculate  the  expected  value  and  the 
variance  of  each  of  the  dependent  variables. 

10.  A  computer  program  has  been  developed  for  the  express  purpose  of 
numerically  evaluating  this  complete  system  of  equations.  Within  this 
program,  computations  are  made  at  successive  times  at  selected  depths  so  that 


the  time  histories  of  the  expected  value  and  the  variance  of  the  dependent 
variables  a  ,  U  ,  and  U  can  be  constructed.  The  program  also  computes 
and  prints  out  the  explicit  contributions  of  each  of  the  input  (independent) 
random  variables  to  the  overall  dispersion  of  the  various  output  quantities 
(dependent  variables) . 


PART  IV:  NUMERICAL  EXAMPLES 
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11.  In  this  section  we  demonstrate  the  application  of  the  above 
methodology  for  conducting  probabilistic  wave  propagation  analyses  for 
bilinear  hysteretic  materials.  Two  numerical  examples  are  presented.  One  is 
referred  to  as  the  "stiff  soil"  problem;  the  other  is  referred  to  as  the 
"soft  soil"  problem.  We  consider  the  former  problem  first.  The  input 
quantities  selected  for  the  first  problem  consisting  of  the  mean  values 

and  the  coefficients  of  variation  for  the  five  independent  varia¬ 

bles  Mo  ’  Mi  *  P  »  PQ  »  and  t  are  given  in  Table  1.  The  coefficients 
of  variation  specified  for  this  problem  are  reasonable  inasmuch  as  they  have 
been  encountered  in  actual  practice.  The  bilinear  material  in  this  example 
problem  is  referred  to  as  "stiff  soil"  because  of  the  high  value  of  the 

loading  modulus  M  .  The  expected  values  E[o  ]  and  E[i)  ]  and  the 
°  max  max 

coefficients  of  variation  v[a  ]  and  v[u  ]  obtained  for  various 

max  max 

depths  from  the  probabilistic  analysis  of  this  problem  are  presented  in 
Tables  1  and  2.  These  tables  also  list  the  individual  contributions  produced 
by  the  constituent  random  variables  X1  to  the  dispersion  of  peak  stress  and 
peak  particle  velocity  at  different  depths.  These  individual  contributions 
are  presented  in  the  last  four  columns  of  Tables  1  and  2  in 


(Table  2).  Note  that  Equation  11  stipulates  that  these  individual  con¬ 
tributions  would  correspond  to  the  coefficient  of  variation  of  peak  stress 
(or  peak  particle  velocity)  if  for  each  variable  X^  the  other  four  constitu¬ 
ent  variables  were  deterministic. 

12.  The  input  quantities  for  the  second  example  problem  are  given  in 
fable  3.  The  bilinear  material  used  for  this  example  problem  is  referred  to  as 
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Table  1.  Individual  Contributions  to  Dispersion  of  Peak  Stress 
Due  to  Uncertainties  in  Input  Variables  (Stiff  Soil;  a  =  0.333) 


|  do 

max 

Si 

DX.  E[ 

0  ] 
maxJ 

xi 

Mi 

VMi 

Z  =  0.0  ft 

Z  =  10.0  ft 

Z  =  20.0  ft 

Z  =  40.0  ft 

M 

63.0  ksi 

0.5 

0.0 

0.01905 

0.03757 

0.07286 

Mi 

252.0  ksi 

0.5 

0.0 

0.00856 

0.01690 

0.03288 

p 

100.0  lb/ft3 

0.1 

0.0 

0.00209 

0.00415 

0.00812 

P 

0.1  ksi 

0.1 

0.1 

0.10008 

0.10010 

0.09994 

T 

0.05  sec 

0.1 

0.0 

0.00420 

0.00831 

0.01627 

E[o  ]  = 

max 

E[o  ]  = 

max 

E[o  ]  = 

max 

E  [o  ]  = 

max 

0.1  ksi 

0.09570  ksi 

0.09168  ksi 

0.08441  ksi 

v[o  ]  = 

max 

V[c  ]  = 
max 

v[o  ]  = 

max 

v[o  ]  = 

max 

0.1 

0.10246 

0.10888 

0.12980 

Table  2.  Individual  Contributions  to  Dispersion  of  Peak  Particle  Velocity 
Due  to  Uncertainties  in  Input  Variables  (Stiff  Soil;  a  =  0.333) 


30 

S. 

max 

l 

1 

3Xi 

e[u  ] 

_ 1 

max 

X. 

i 

VMi 

Z  =  0.0  ft 

Z  =  10.0  ft 

Z  =  20.0  ft 

Z  =  40.0  ft 

M 

o 

63.0  ksi  0.5 

0.23618 

0.21932 

0.20262 

0.17008 

Mi 

252.0  ksi  0.5 

0.0 

0.00779 

0.01546 

0.03033 

p 

100.0  lb/ft3  0.1 

0.04544 

0.04760 

0.04970 

0.05372 

P 

o 

0.1  ksi  0.1 

0.09074 

0.09123 

0.09164 

0.09229 

T 

0.05  sec  0.1 

0.0 

0.00384 

0.00761 

0.01502 

e[u  ]  = 

max 

e[u  ]  = 

max 

e[u  ]  = 

max 

e[u  ]  = 

max 

2.99  fps 

2.85  fps 

2.72  fps 

2.48  fps 

v[u  ]  = 

max 

v[u  ]  = 

max 

v[u  ]  = 

max 

v[u  ]  = 

max 

0.25706 

0.24266 

0.22898 

0.20447 

"soft 

soil"  because  of  the 

low  value 

of  the  loading  modulus  M 

0 

.  The  pro- 

babilistic  output  quantities  for  this 

problem  are 

presented  in 

Tables  3  and 

13.  The  information  presented  in  Tables  1  through  4  reveals  several 
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Table  3.  Individual  Contributions  to  Dispersion  of  Peak  Stress 
Due  to  Pncer tainties  in  Input  Variables  (Soft  Soil;  a  =  0.714) 


3a 

max 

Si 

3Xi 

E[ 

0  1 
maxJ 

X. 

1 

Mi 

Si/lJi 

Z  = 

0.0  ft 

Z  = 

10.0  ft 

Z  =  20.0  ft 

Z  =  40.0  ft 

M 

o 

7.0  ksi 

0.5 

0.0 

0. 

.04475 

0.08463 

0.14888 

Mi 

252.0  ksi 

0.5 

0.0 

0. 

,00246 

0.00451 

0.00751 

p 

100.0  lb/ft3 

0.1 

0.0 

0. 

00819 

0.01559 

0.02777 

P 

o 

0.1  ksi 

0.1 

0.1 

0. 

,10153 

0.10263 

0.10377 

T 

0.05  sec 

0.1 

0.0 

0. 

,01640 

0.03122 

0.05555 

E[o  ]  = 
max 

E[a  ]  = 

max 

E[o  ]  = 

max 

E[a  ]  = 

max 

0.1 

ksi 

0.0834  ksi 

0.0706  ksi 

0.0525  ksi 

v[a  ]  = 

max 

v[o  ]  = 
max 

v[o  ]  = 

max 

v[o  ]  = 
max 

0.1 

0.11249 

0.13760 

0.19195 

Table  4. 

Individual  Contributions 

to  Dispersion 

of  Peak  Particle  Velocity 

Due 

to  Uncertainties  in  Input 

Variables  (Soft  Soil;  a  = 

0.714) 

aii 

max 

Si 

E[U  ] 

L  maxJ 

X. 

1 

P 

i 

VMi 

Z  =  0.0  ft 

Z  = 

10.0  ft 

Z  =  20.0  ft 

Z  =  40.0  ft 

M 

o 

7.0 

ksi 

0.5 

0.23618 

0. 

19938 

0.16526 

0.10806 

Mi 

252.0 

ksi 

0.5 

0.0 

0. 

,00211 

0.00391 

0.00657 

P 

100.0 

lb/ft3 

0.1 

0.04544 

0. 

,05403 

0.06186 

0.07470 

P 

o 

0.1 

ksi 

0.1 

0.09074 

0. 

09295 

0.09479 

0.09723 

T 

0.05 

►  sec 

0.1 

0.0 

0. 

,01497 

0.02878 

0.05195 

E[U  ]  = 
max 

E[U  ]  = 

max 

e[u  ]  = 

max 

e[u  ]  = 

max 

8.98  fps 

7.43  fps 

6.23  fps 

4.57  fps 

v[u  ]  = 

max 

v[u  ]  = 

max 

v[u  ]  = 

max 

v[u  ]  = 

max 

0.25706 

0.22703 

0.20240 

0.17162 

interesting  trends  regarding  the  dispersion  of  peak  stress  and  peak  particle 


velocity: 


a.  For  the  two  soil  materials  considered  in  this  paper,  the 

coefficient  of  variation  of  peak  stress  \'[a  1  increases 

with  depth.  On  the  other  hand,  the  coefficient  of  variation 
of  peak  particle  velocity  v[u  x]  decreases  with  depth.  The 
rates  of  change  of  the  coefficients  of  variation  with  depth, 
however,  are  more  pronounced  for  the  soft  soil  case. 

jj.  As  expected,  at  the  ground  surface  only  the  variability  in 
peak  pressure  Pq  contributes  to  the  dispersion  of  peak 
stress.  With  increasing  depth,  however,  the  variabilities  in 
the  other  input  variables  also  contribute  to  the  dispersion  of 
peak  stress.  Moreover,  the  contribution  of  the  variability  in 
is  (for  all  practical  purposes)  the  same  for  all  depths 
whereas  the  influence  of  the  other  input  variables  becomes 
more  significant  with  increasing  depth.  These  influences  are 
stronger  for  the  soft  soil  case.  For  both  soils,  the  varia¬ 
bilities  in  peak  pressure  P^  and  loading  modulus  >1q 
contribute  most  to  the  dispersion  of  peak  stress.  The  combined 
variability  in  unloading  modulus  ,  density  o  ,  and  decay 
constant  t  has  a  very  small  influence  on  the  dispersion  of 
peak  stress. 

c_.  The  variability  in  loading  modulus  M  contributes  most  to 

dispersion  of  peak  particle  velocity.  Its  influence,  however, 
decreases  with  increasing  depth.  The  variability  in  unloading 
modulus  has  a  negligible  effect  on  the  dispersion  of  peak 

particle  velocity  but  its  influence  increases  with  increasing 
depth.  The  combined  variability  in  density  c  ,  peak  pressure 
Pq  ,  and  decay  constant  I  has  a  small  effect  on  the 

dispersion  of  peak  particle  velocity  but  its  influence 
increases  with  increasing  depth  and  is  more  pronounced  for  the 
soft  soil  case. 

14.  It  is  worth  pointing  out  that  for  the  two  example  problems 

considered,  the  contributions  of  the  second  partial  derivatives  in  Equation 

10  to  the  expected  values  of  the  dependent  variables  a  ,  U  ,  U  ,  and  t  are 

a 
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not  very  significant.  In  the  case  of  stress,  the  contribution  is  negligible 
For  the  other  three  dependent  variables  the  contributions  on  the  average  are 
about  10  percent  of  the  expected  values  and  as  high  as  25  percent  for  the 
particle  velocity  at  late  times.  The  contributions  are  more  pronounced  for 
the  soft  soil  problem. 

15.  To  provide  a  feeling  for  the  effects  of  input  variability  on  the 
dispersion  of  particle  displacement  for  the  two  problems  considered,  the 
results  of  the  probabilistic  calculations  for  surface  displacement  at 
t  =  0.15  second  are  presented  in  Table  5.  This  tabulation  shows  that  the 
variability  in  loading  modulus  Mq  obviously  contributes  most  to  dispersion 
of  particle  displacement  for  both  of  the  example  problems.  The  combined 
effect  of  the  variabilities  in  unloading  modulus  ,  density  p  ,  peak 

pressure  Pq  ,  and  decay  constant  r  on  the  dispersion  of  particle  displace 
ment  is  small. 


Table  5.  Individual  Contributions  to  Dispersion  of  Particle  Displacement 
Due  to  Uncertainties  in  Input  Variables  (Z  =  0.0  ft;  t  =  0.15  sec) 


Soft  Soil,  ct=  0.714 _  _ Stiff  Soil,  a  =  0.333 


3U 

s. 

1 

3U 

si 

Xi  Vyi 

E[U]  p.  S./p. 

3X. 

l 

E[U] 

M 

o 

7.0  ksi 

0.5 

0.29385 

63.0  ksi 

0.5 

0.29738 

M 

252.0  ksi 

0.5 

0.05906 

252.0  ksi 

0.5 

0.05999 

P 

100.0  lb/ft3 

0.1 

0.04459 

100.0  lb/ft3 

0.1 

0.04595 

P 

o 

0.1  ksi 

0.1 

0.08905 

0.1  ksi 

0.1 

0.09175 

T 

0.05  sec 

0.1 

0.02215 

0.05  sec 

0.1 

0.04232 

E[U]  =  E[U]  = 

10.92  in.  2.73  in. 


V[U]  = 
0.31661 


V[U]  = 
0.32303 


16.  Typical  time  histories  (wave  forms)  of  the  expected  value  and  the 
coefficient  of  variation  of  stress,  particle  velocity,  and  particle  dis¬ 
placement  at  Z  =  40.0  feet*  are  portrayed  in  Figures  1  through  3,  respec¬ 
tively.  The  arrival  times  for  these  graphs  correspond  to  Eft^]  =  0.0257 

second  for  the  stiff  soil  and  Eft  1  =  0.0770  second  for  the  soft  soil.  The 

a 

coefficient  of  variation  of  the  arrival  time  is  a  function  only  of  the 
variability  in  the  independent  variables  and  p  and  is  equal  to 

V[t  ]  =  0.2416  for  all  depths  and  for  both  problems  considered.  The 
individual  contributions  to  the  dispersion  of  t  due  to  uncertainties  in 

and  p  are,  respectively,  0.2373  and  0.0456.  Therefore,  the  variability 
in  Mq  contributes  most  to  the  dispersion  of  the  arrival  time.  It  is 
interesting  to  note  from  Figures  2  and  3  that  the  coefficients  of  variation 
of  particle  velocity  and  particle  displacement,  v[ll]  and  v[u],  respec¬ 
tively,  also  increase  with  increasing  time.  In  the  case  of  stress  (Figure  1) 
the  coefficient  of  variation  v[o]  initially  decreases  with  time  and  then 
increases  at  late  times  (the  soft  soil  history  eventually  turns  upward  at 
t  >  0.2  second).  The  increase  in  v[u]  with  time  is  important  since  the 
expected  value  of  particle  displacement  E[u]  also  increases  with  time. 

17.  Information  of  the  type  presented  in  Tables  1  through  5  and 
Figures  1  through  3  can  be  used  to  perform  probabilistic  analysis  of  airblast 
induced  ground  shock  and  to  construct  reliability  based  safety  factors  for 
blast-resistant  design  purposes.  Assuming  that  the  probability  distribution 
functions  for  each  of  the  dependent  random  variables  are  roughly  bell-shaped, 
the  normal  distribution  function  can  be  used  to  conduct  the  desired  proba¬ 
bility  analysis.  For  example,  applying  the  rule-of-thumb  that  the  proba¬ 
bility  that  a  variable  lies  within  two  standard  deviation  bounds  of  its  mean 
is  approximately  95  percent  to  the  surface  displacement  predictions  in 


Table  5,  it  can  be  stated  P[4.006  in.  _<  U  £  17.834  in.]  =  0.95 
for  the  soft  soil  and  P[0.966  in.  _<  U  4.493  in.]  =  0.95  for  the  stiff 
soil.  In  the  absence  of  the  exact  shape  of  the  probability  distribution 
function  one  could  also  use  the  Chebyshev  inequality  (Reference  2)  for 
probability  analysis. 
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PART  V:  CONCLUSIONS 


18.  The  partial  derivative  method  has  been  used  to  convert  a  determi¬ 
nistic  analytic  solution  for  stress  wave  propagation  in  homogeneous  bilinear 
hysteretic  materials  subjected  to  an  exponentially  decaying  surface  airblast 
pulse  into  a  numerical  probabilistic  solution.  The  probabilistic  solution 
has  been  coded  and  used  to  conduct  sensitivity  analyses  to  rank  the  relative 
effects  of  input  variabilities  and  uncertainties  on  the  dispersion  of  output 
quantities  (stress,  particle  velocity,  and  particle  displacement).  The 
analyses  indicate  that  variability  in  the  material  loading  modulus  contributes 
most  to  the  dispersion  of  particle  velocity  and  particle  displacement  but  its 
influence  decreases  with  increasing  depth.  The  variabilities  in  loading 
modulus  and  peak  airblast  pressure  contribute  most  to  the  dispersion  of  peak 
stress;  the  influence  of  the  former  increasing  with  increasing  depth  and  that 
of  the  latter  remaining  essentially  constant  with  depth.  The  analyses  also 
indicate  that  the  coefficient  of  variation  of  particle  velocity  increases 
with  time  whereas  the  coefficient  of  variation  of  stress  initially  decreases 
with  time  and  then  increases  at  later  times.  At  any  given  depth,  of  course, 
the  expected  values  of  stress  and  particle  velocity  decrease  with  time.  The 
coefficient  of  variation  of  particle  displacement  increases  with  time,  which 
is  of  particular  interest  since  its  expected  value  also  increases  with  time. 
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APPENDIX  A 
NOTATION 


Cq  Loading  wave  speed 

Unloading  wave  speed 

EL  ]  Expectation  of  a  random  variable 

k  Number  of  standard  deviations  above  and  below  mean  at  which  the 
dependent  random  variable  is  evaluated 

M  Initial  loading  modulus 

Unloading  modulus 

P(t)  Applied  surface  pressure-time  history 

Po  Peak  applied  pressure 

P[  ]  Probability  of  a  random  variable 

S.  Standard  deviation  of  X. 

1  1 

t  T  ime 

t  Arrival  time  =  Z/C 

a  o 

U  Particle  displacement 

U  Particle  velocity 

U  Peak  particle  velocity 

max 

V[  ]  Coefficient  of  variation  of  a  random  variable 
Var[  ]  Variance  of  a  random  variable 

X^  Functionally  independent  random  variable 

Y  Functionally  dependent  random  variable 

Z  Depth 

or  - 1 

a  Hysteretic  parameter  =  -  —  -  - - 

+  1 


Exponential  time  constant  (time  at  which  the  applied  surface 
pressure  has  decayed  to  0.3678  P  ) 

u .  Mean  of  X . 

1  1 

p  Mass  density 

(v  Stress 

a  Peak  stress 

max 
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